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Abstract. The near-infrared behavior of the red giant branch (RGB hereafter) as a function of abundance is 
examined with an unprecedented large sample of 27 Galactic globular clusters with Two Micron All Sky Survey 
photometry. We propose a new simplified analyses, involving the zero point of the RGB slope fit, and derive 
calibrations for the RGB slope, zero point, and tip. The weak metallicity sensitivity of the zero point leads to a 
"fan"-like diagram for obtaining the abundance distributions in resolved stellar systems, and reddening estimates. 
Finally, we apply the new calibrations to the recently discovered Galactic globular cluster 2MASS GC01, to 
derive [Fe/H]H96 = —1.19 ±0.38 mag. The uncertainty is dominated by the severe foreground contamination. We 
estimate an extinction of Ay = 21.07 ± 2.20 mag toward the cluster. 
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1. Introduction 

The red giant branch (RGB hereafter) stars are among 
the brightest red stars in stellar systems, older than a few 
Gyrs. These stars appear in almost all galaxies, including 
II Zw 40 (Ostlin |2000 ), considered until recently as the 
best candidate for a primeval galaxy. Therefore the red 
giants are a promising tool for probing the parameters of 
old populations and the history of star formation in any 
galaxy. 

Galactic globular clusters, with their single age and 
metallicity, are the ideal sites for calibrati ng the RGB pa- 
rameters. Since Da Costa & Armandroff ( 1990 ) provided 
the first reliable calibration of the position of the RGB ver- 
sus metallicity, there has been a significant advancement, 
both because of the improvement of the astronomical in- 
strumentation and the development of the corresponding 
theory. 

The infrared waveband is particularly compelling for 
such studies, in comparison with the optical one, because 
of the potential to probe the stellar populations of sys- 
tems with high foreground and/or intrinsic extinction. 
The relatively small size and field of view of the IR ar- 
rays have made it more difficult to carry out photometry 



of large areas, and to compile uniform samples, necessary 
to calibrate the RGB parameters reliably, in comparison 
with the optical region. The Two Micron All Sky Survey^ 
(2MASS hereafter) offered for the first time such an op- 
portunity. 

Previous calibrations suffered a number of drawbacks. 



Kuchinski, Frogel, & Terndrup (1995) and Kuchinski & 
Frogel ( 1995| ) studied only metal rich globular clusters, 
with [Fe/H]> —1. Later on, Ivanov et al. ( 2000a ) added 
to their sample three metal poor globular clusters. Ferraro 
et al. ( |2000| ) used exceptional quality data but in a pho- 
tometric system, based on unpublished standards by Ian 
Glass (South African Astronomical Observatory) with no 
available transformations to any of the other systems. The 
similarity of their final results to those of Ivanov et al., who 
used the CIT system, leads to a conclusion that the two 
photometric systems are not radically different. 

To expand the basis of the RGB studies, we used high 
quality uniform photometry of Galactic globular clusters 
from the 2MASS point source catalog, and assembled a 
sample of RGBs of clusters with well known distances and 
reddening. We calibrated the behavior of the RGB slope, 
zero point, and tip with metallicity in the 2MASS photo- 
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metric system. For the first time we offer such calibration 
in a well defined photometric system with all-sky cover- 
age. We present a "fan" diagram, suitable for abundance 
distribution analysis. Our results complement the recent 
work of Cho & Lee (2001) and Grocholski & Sarajedini 



(2001). who explored the properties of the RGB bump. 



2. 2MASS Calibration of the RGB Parameters 

2.1. Sample 

Infrared photometry of about 80 of the 147 globular clus- 



ter, listed in Harris (1996, revision June 22, 1999) is cur- 
rently available from the First and the Second 2MASS 
incremental releases. However, many of the observed clus- 
ters are rendered unsuitable for this project because they 
suffer from one or more of the following drawbacks: 

1. Severe foreground contamination. This refers mostly 
to globular clusters in the general direction of the 
Galactic center. 

2. The photometry of some distant halo clusters or highly 
reddened clusters is not sufficiently deep. 

3. Compact globular clusters have unresolved cores, lim- 
iting the number statistics, and leading to crowding 
problems. 

After a visual inspection of the luminosity functions 
and color-magnitude diagrams (CMD hereafter) of all 
available globular clusters, we selected a subset of 27 ob- 
jects (Table |l|) for which reliable estimates of at least one 
RGB parameters could be obtained. To calculate the ab- 
solute magnitudes, and the intrinsic colors of stars, we 
adopted the reddening and distance estimates from the 
compilation of Harris (1996), and the reddening law of 
Rieke & Lebofsky (1985; used throughout this paper). 



A special care was taken t o ass ure the uniformity of 
the metallicity data. Harris (1996) based his co mpila - 
tion o n the system established by Zin n fc Wsst ( 1984 ), 
Zinn ( |1985[), a nd Armandroff & Zinn ( |1988|) . Carretta & 
Gratton ( 1997 ) suggested that this scale may overestimate 
the abundance of metal-rich globular clusters, and devel- 
oped a new one. We list their measurements in Table ^ as 
well. Values in brackets are estimates based on the trans- 
formation they derive, to the Zinn & West scale (see Eq. 
7 of Carretta & Gr atton ). 

Ferraro et al. ( |1999|) argued that the total amount 
of heavy elements is a better metallicity measurement, 
because it accounts naturally for the opacity variations, 
which depend on the total metallicity, not just the iron 
abundance (Salaris, Chieffi, & Straniero 1993). They de- 
veloped a news scale, designed to measure the that. For 
the clusters, absent in their sample, we adopted their pre- 
scription for ev-element enhancement (Sect. 3.4 in Ferraro 
et al. ; estimates bracketed in Table [l]) . 

Since the metallicities originate from different sources, 
and the corresponding uncertainties are often not quoted, 
we adopted an uniform error of 0.20 dex for all measure- 
ments. This is a typical value of the accuracy of the abun- 
dance estimates. 



Table 1. Globular cluster sample. The Columns give: 
cluster ID, color excess Eb-v, distance modul us, tid al ra- 
dius Rt (all from the compilation of Harris 1996| ), and 
abundances in various scales (see the references). 
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References for [Fe/H|: 



H96 - Harris ( |1996j ; Zinn scal e); 
CG97 - Carretta & Gratto n jl997 
F99 - Ferraro et al. (|l999|) 



The first step before estimating the RGB parameters 
was to eliminate statistically the foreground contamina- 
tion. As many stars were removed randomly from the 
CMD of the globular clusters, as the number of stars, 
present in a nearby field with the same area as the glob- 
ular cluster field. The statistical removal is not reliable in 
case of severely contaminated globular clusters, where the 
number of foreground stars was comparable to the number 
of cluster members. A subset of CMDs for four clusters is 
shown in Figure [l] and |[ In addition, we removed the stars 
within 15 arcsec from the clusters centers, to minimize the 
crowding effects. 

2.2. RGB Slope and Zero Point 



Da Costa & Armandroff (1990) prepared a set of stan- 
dard RGBs on the V versus V-I CMD, for metallicity es- 
timates of globular clusters in the optical. The technique 
has been expanded toward the near infrared wavebands, 
where the giant branch is linear, and therefore, more ob- 
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Fig. 1. Color-magnitude diagrams for a subset of metal 
poor globular clusters. The left panels show the raw data. 
The right panels show the reddening corrected CMDs, 
converted to absolute magnitudes. The foreground con- 
tamination and stars with (J — Ks)o < 0-4 mag have 
been removed. 
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Fig. 2. Same as Figure [j] for metal rich clusters. 



vious and easier to define compared with the optical or 
optical-infrared CMDs. The slope is related to the effec- 
tive temperature of the stars along the RGB, and T e ff in 
turn depends on the opacity, and the heavy element abun- 
dance. The slope is a reddening and distance-free metal- 



licity indicator. It is significantly less demanding in terms 
of observing time and telescope collecting area than the 
spectroscopic methods, but unlike them it can be applied 
reliably only to uniform groups of stars. 

We followed the procedure to determine the RGB 



slope, described in Ivanov et al. (2000a), fitting the RGB 
on the [K a , (J — K s )] CMD with a linear equation: J — 
K s = a x K s + b. Only the stars above the horizontal 
branch were included. We applied a least square method, 
taking into account the uncertainties along both axis. Two 
iterations were carried, and all the stars outside 5<r from 
the first fit were removed from the calculations. 

For the purpose of fitting the RGB of an individual 
cluster, the zero point is equivalent to previous calibra- 
tions of the RGB c olors a t fixed if-band magn itude levels 
(e.g. Ivanov et al. 2000a ; Ferraro et al. 200C), although 
more straightforward. To derive the zero-point relation to 
the abundance, we used only the clusters which suffer min- 
imal extinction (Eb-v <0.3 mag). 

Notably, the fit to the zero point shows a very small 
variation with metallicity: 0.11 mag for [Fe/H]#g6 varying 
from -2 to -0.5. For comparison, (J — K$)o at Kg = —5.5 
mag varies by 0.32 mag for the same metallicity range 



(Ferraro et al. 2000). This result offers the possibility for 
a more reliable reddening estimates than in case of RGB 
color calibrations. 

The behavior with metallicity of the RGB slope and 
zero point are demonstrated in Figure Fitting coeffi- 
cients for individual cluster RGBs are listed in Table [| 



2.3. RGB Tip 

The tip of the RGB is a well known distance indicator. 
Mould & Kristian fll986|), Lee Ql993|) and Lee et al. |1993 



pioneered the method in the optical, measuring the cor- 
responding jump in the luminosity function in external 
galaxies. Later it was applied to number of resolved stellar 
systems, among the most challengin g of w hich were some 
Virgo cluster mem bers (Harris et al. 1998 ) and NGC 3379 
(Sakai et al. |l997| ). 

This is a statistically demanding technique, requiring 
50-100 stars per bin. The limited size of infrared arrays 
explains the difficulty to apply it in the near infrared. 
Although it does cover a large area, the 2MASS photom- 
etry can not alleviate the intrinsic problem of the small 
number of giants in globular clusters. Thus, we assume 
that the brightest cluster member represents the RGB tip. 

In most of the cases the tip is obvious, but sometimes 
additional criteria had to be applied to determine the 
brightest stars. We took advantage of the linearity of the 
RGB in the near infrared, and excluded from the consid- 
erations bright stars that deviated from the color of the 
RGB at a given magnitude level, predicted by the RGB 
slope fit by more than 0.5 mag. We also excluded some ex- 
tremely bright stars, with luminosity higher than the rest 
of the RGB by 2-3 mag. They were obvious foreground 
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Table 2. Linear fits to the RGB in the form of: J — K a = 
axK s + b. The radii R max of the region around the cluster 
center used in the fit, and the root mean square of the fit 
are given. 
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Points for NGC 6273 and 6656 are omitted because 
of the large reddening toward this cluster. 
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Fig. 3. Relation of the [Fe/H] in various metallicity scales 
versus RGB slope on [Mk 3 , (J — K s )q] diagram. The bars 
indicate la u ncert ainties. The dashed line is the fit of 
Ferraro et al. ( |2000[) drawn without transformation of col- 
ors. The zero point errorbars show only the statistical er- 
rors from the fits to the RGB, omitting the 0.2 mag added 
in quadrature to account for the uncertain distance moduli 
and reddening corrections. 



contamination. Finally, the red variables from Clement et 
al. ( pOOll ) were excluded . 

The formal uncertainties of the stellar magnitudes 
given in 2MASS were discarded, since they do not repre- 
sent well the uncertainty in the tip magnitude. Instead, we 
adopted the difference in the magnitudes of the two bright- 
est stars, accounting for the possibility that the brightest 
star may be a non-member. This led to typical error values 
of 0.2-0.4 mag, much larger than the 2MASS errors. 

The RGB tip magnitudes for 20 globular clusters are 
given in Table ||. The behavior of the RGB tip with metal 
abundance is shown in Figure ^. 

2.4. Results 

The behavior of the derived RGB parameters with metal- 
licity was fitted with linear equations, taking into account 
the errors along both axes. The coefficients are given in 
Table |. Figure | shows a "fan" -like grid of RGBs in a red- 
dening and distance corrected CMD, for the three metal- 



licity scales discussed in Section 2.1. It demonstrates that 



the RGB behavior can be reduced to a simple rotation 
around a nearly-fixed point. This comes to no surprise, 
since the zero point of the RGB fit is almost independent 
of the abundance (Section |2.2|) . 
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Fig. 4. Relation of the [Fe/H] in various metallicity scales 
versus RGB tip absolute magnitude. The bars indicate ltr 
uncertainties, with the Y-axis errors including the uncer- 
tainties due to the reddening and distance. The dashed 
line is the fit of Ferraro et al. (|2000D drawn without trans- 
formation of colors. 
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Table 3. Estimated RGB tip absolute magnitudes. 



NGC 



J 



H 



K s 



Table 4. Polynomial fits to the metallicity behavior of the 
RGB slope, zero point, bump, and tip: Y=ao+aixX. 



104 
288 
1851 
1904 
5024 
5139 
6144 
6171 
6205 
6284 
6441 
6624 
6637 
6656 
6779 
6809 
6864 
7078 
7089 
7099 



616(0.171) 
578(0.621) 
182(0.083) 
055(0.189) 
093(0.049) 
124(0.106) 
062(0.294) 
326(0.605) 
156(0.050) 
267(0.277) 
521(0.246) 
459(0.019) 
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146(0.023) 
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-6.073(0.133) 
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-6.862(0.209) 
-6.672(0.656) 
-6.306(0.040) 
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Note: The difference between the two brightest 



stars is given in brackets. 



The grid allows to obtain the metallicities of individual 
sdtars in resolved systems in the infrared, and to obtain 
the metallicity distributions as it has been done before in 
the optical (e.g. Saviane et al. 2000| ). The linear represen- 
tation of the RGB in absolute Ks magnitude versus the 
intrinsic (J — K$)o color is: 



M Ks = RGB ZP + RGB S i x ( J - K s ) 



(1) 



Both the slope RGB si and the zero point RGBzp were 
calibrated as functions of the abundance: 



RGBsi = a 
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[Fe/H] 
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(2) 
(3) 



After substituting Eq. 2 and 3 in 1, and solving for the 
stellar metallicity: 

[Fe/H] = [M Ks -a sl x(J-K s ) -a$ p } / [af l x (J-K s )o+a 

This equation should be used with caution. Small pho- 
tometric errors are crucial for accurate [Fe/H] estimates. 
The upper part of the infrared CMD is better suited for 
determining individual metallicities because near the root 
of the RGB the abundance relation degenerates. The ex- 
act faint limit is determined by the quality of the pho- 
tometry, and by the presence of the horizontal branch at 
Mk s ~ —1.5 mag. 

3. 2MASS GC01 
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-6.712(0.210) 
-6.631(0.193) 
-7.032(0.212) 
-7.109(0.243) 
-7.003(0.224) 



-0.051(0 
-0.058(0 
-0.060(0 
-0.070(0 
-0.082(0 
-0.084(0 
-0.323(0 
-0.387(0 
-0.399(0 
-0.486(0 
-0.594(0 
-0.620(0 
-0.615(0 
-0.739(0 
-0.768(0 



006 
007 
008 
104 
120 
122 
121 
151 
158 
121 
153 
162 
134 
171 
181 



002 
002 
002 
007 
007 
007 
025 
026 
026 
025 
024 
024 
032 
034 
034 



22 
22 
22 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 



Note: Errors are given in brackets. Metallicity 
notation is the same as in Table [l| 




Hurt ct al. (1999, 2000| ) reported the serendipitous discov- 
ery of 2MASS GC01 (hereafter GC01) in 2MASS data. 



0.4 0.6 0.8 1 1.2 

(J-K s ) 

f P ](4) 

Fig. 5. "Fan Diagram" - a grid of RGBs on the reddening 
and distance corrected color-magnitude diagram, for the 
three metallicity scales, as indicated. 



It is a heavily obscured globular cluster, laying in the 
Milky Way disk, in the general direction of the Galactic 
center. Ivanov, Borissova, & Vanzi ( 2000b| ) determined 
(m - M) = 12.4 - 14.0, and Ay = 20.9 - 18.8, assuming 
[Fe/H] = — 0.5 and —2.0, respectively. The main source of 
uncertainty in these estimates was the unknown metallic- 
ity of the cluster, although the location of GC01 suggested 
that it might be a metal rich object. 
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The first step toward a metallicity estimate of GC01 
was to remove the foreground stars contamination. Unlike 
the clusters we used to derive the RGB parameters calibra- 
tions, the contamination here is severe, reaching ~ 30% in 
the RGB region. We performed 2000 foreground subtrac- 
tions, and estimated the RGB slope for each realization 
separately. This method yields distributions of the RGB 
slope and tip, and the widths of these distributions mea- 
sure the respective uncertainties. 

To carry out this procedure we defined the CMD area, 
encompassing the RGB: 8.0 < K$ < 14.0 mag, and 
2.5 < J - K s < 6.0 mag. Then we divided it into 0.2 
mag square bins. Experiments with different bin sizes in- 
dicated that any value between 0.2 and 1.0 mag leads to 
the same conclusions. 

Next, we constructed CMDs for the cluster field, and 
for a surrounding field with an equal area. To minimize the 
crowding effects we omitted the stars within 15 arcsec from 
the cluster center. The outer limit of the cluster field was 
const rained by th e cluster diameter (3.3±0.2 arcmin, Hurt 
et al. 1999 , 2000 ). We carried out our calculations to two 
values of the outer radii: 1.0 and 1.5 arcmin. Smaller values 
limit the number statistics, and larger ones increase the 
fraction of the foreground contamination. The foreground 
field was defined as a circular annulus with 10.49 — 10.54 
or 10.43 — 10.54 arcmin size, respectively for 1.0 and 1.5 
armin cluster fields. 

Finally we counted the stars in each bin, and sub- 
tracted randomly from the "cluster" bins as many stars 
as were present in the "field" bins. If the latter bin had 
more stars than the former one, we subtracted stars from 
the nearby bins, again in a random way. 

We carried out a linear fit on the RGB stars in the 
foreground subtracted CMD in th e sa me manner as for 
the calibration clusters (see Section 2/2). For stars without 
error measurements in 2MASS we adopted a = 0.20 mag. 
Using the linearity of the RGB, we imposed a faint limit 
of the stars, included in the fit, just above the horizontal 
br anch lev el. The luminosity function of GC01 (Ivanov et 
al. 2000b , Figure 5) indicates that the horizontal branch 
is at -?Cs=13.0-13.2 mag. To minimize the uncertainties of 
the RGB slope, we also imposed color limits on the stars 
we used in the fit. The red one was set to J— K$=5.5 mag, 
and has no effect on the slopes because of the negligible 
number of stars to the red of the RGB. The results are 
somewhat more sensitive to the blue limit. We choose to 
impose J — Ks—3-50 and 3.80 mag, because of the well- 
defined limit of the RGB at this colors (Figure ||). 

A summary of the results for the RGB slope and zero 
point GC01 is presented in Figure 0. The RGB tip is 
omitted because it is sensitive only to the adopted cluster 
radius. Clearly, the effects from the assumed parameters 
are smaller or comparable with the uncertainties origi- 
nating due to the foreground contamination. We calcu- 
lated the error-weighted averages for the realizations with 




Fig. 6. 2MASS color-magnitude diagram of GC01. Solid 
dots are stars from 0.25 to 1 arcmin from the cluster center 
- cluster+foreground. X's indicate the stars within 10.49 
to 10.54 arcmin from the center - pure foreground, la 
errors are indicated. The solid line represents the average 
RGB fit, described in Section |[ Clearly, it is dominated 
by the stars with minimal uncertainties. Dotted lines are 
the adopted limits. 



1. < Slopc>= -0.0965 ±0.0171 

2. < ZeroPoint > = 5.40 ± 0.18 

3. < Tip > = 8.87 ±0.13 

Using the equations, given in Table [|, we obtain: 
[Fe/H] ff96 = -1.19±0.38, [Fe/H] CG97 = -1.06±0.34, and 
[M/H] F99 = -0.88±0.32. The errors given here include 
both the uncertainties in the slope, and in the calibra- 
tions. GC01 appears to resemble closely NGC 104. 

The estimated absolute magnitude for the RGB tip, for 
a cluster with such abundance is Mk s = — 6.30±0.35 mag, 
using our new calibration for the Zinn metallicity scale. 
This leads to a distance modulus of (m — M) p + A y, = 
15.17±0.38 mag consistent with Ivanov et al. ( [20001 ). We 
refrain from further considerations based on the RGB tip 
because of the poor statistics at the brighter end of the 
RGB. Instead, we will adopt (m - M) + A Ks = 15.0 ± 
0.4 mag, a result of interpolation between the values for 
[Fe/H]=-1.0 and -2.0, in Table 2 of Ivanov et al. (§000tj. 

We can also verify if the cluster reddening is consistent 
with the previous estimates. First, we subtract the RGB 
slope equations, written for the GC01 in apparent, and in 
absolute magnitudes. Respectively: 



T.~t i in 



!,m = 12.5 — 13.0, to avoid the possible influence of the 
cluster horizontal brunch. The determined RGB parame- 
ters for GC01 are: 



(mj - m Ks ) = RGB^p + RGB S1 x m Ks 



(Mj - M Ks ) = RGB|p s + RGB S i x M Ks 



(5) 
(6) 
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Fig. 7. Behavior of the RGB slope (top) and zero point 
(bottom) for different parameters of the foreground sub- 
traction and fits. The horizontal axis is the lower limit of 
the stars used in the RGB fit. Triangles indicate cluster 
radius of 1.0 arcmin, and circles indicate 1.5 arcmin. Open 
symbols indicate a blue limit of J — i^5=3.25 mag, and 
solid symbols of J — Kg=3.00 mag. Only some la errors of 
individual points are shown to avoid crowding. The solid 
lines indicate the weighted averages, and the dashed lines 
indicate their la errors. The values to the right of the ver- 
tical dotted line were discarded because they were affected 
by the horizontal branch. 

Here the notation is the same as in Table |^, and RGB|p s = 
0.36 ±0.20 mag, for the [Fe/H]# 96 given above. The result 
of the subtraction is: 



diagram, suitable for analyses of the metallicity spread in 
resolved stellar systems. 

The derived calibrations were applied to estimate 
the metal abundance of the recently discovered globu- 
lar cluster GC01. It is a particularly challenging object 
because of the severe foreground contamination. We re- 
moved it with a random procedure, and used the RGB 
slope of the remaining pure cluster population to derive 
[Fe/H]#96=— 1.19±0.38. The uncertainty is dominated by 
the foreground contamination, and albeit large, allows to 
exclude extremely high abundances, expected from the 
cluster location. GC01 is likely to be moderately metal 
poor. The RGB tip and zero point yield distance modulus 
and extinction, consistent with our previous estimates. 
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